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ABSTRACT 



Aims. We explore the color dependence of the radial profile of satellite galaxies around isolated parent galaxies. 

Methods. Samples of potential satellites selected from large galaxy redshift surveys are significantly contaminated by interlopers - 

objects not bound to the parent galaxy. We use the Sloan Digital Sky Survey to estimate the interloper fraction in samples of candidate 

satellite galaxies. 

Results. We show that samples of red and blue satellites have different interloper populations: a larger fraction of blue galaxies are 

likely to be interlopers compared to red galaxies. Both with and without interloper subtraction, the radial profile of blue satellites 

is significantly shallower than that of red satellites. In addition, while red and blue primaries have different interloper fractions, the 

slope of the corrected radial profiles are consistent after interloper correction. We discuss the implications of these results for galaxy 

formation models. 

Key words. cosmology:theory - dark matter - galaxies: formation - galaxies: fundamental parameters - galaxies: structure 



1. Introduction 

In cold dark matter (CDM) models of the universe, large num- 
bers of dark matter (DM) subhalos lie within the virial radius 
of larger dark matter halos. In some halos and subhalos, bary- 
onic material has cooled and formed stars, resulting in a central 
galaxy and satellite galaxies. The spatial distribution of satel- 
lite galaxies in galaxy-sized halos, then, reflects the evolution 
of satellite galaxies and the mass accretion history of their par- 
ent halos. For example, dark matter simulations suggest that the 
angular distribution of subhalos follows the shape of the DM 
halo, which is indicative of infall of subhalos along filaments 
(e.g.. lZentner et al]|2005l : lLibeskind et al.l l2005), while observa- 
tions suggest that satellites lie along the major axis of the light 
distribution for early-type galaxies (e.g.. Sales & Lambasll2004t 



lBrainerdll2005l: lYang et al.lbooel: lAzzaro et al.ll2007l) 

In dark matter simulations, the radial distribution of subhalos 
is biased with respect to the density profile of DM halos; at small 
separations from the center of the halo (within ~ 20 - 50% of the 
virial radius), the distribution of DM subhalos has a lower con- 
centr ation, bu t it follows the dark mat ter density profile at larger 
radii dGhigna et al. 1998; Colin e t al.' 1999; Ghigna et al. 2000; 
Springel et alJl2001l;lDe Lucia et al. 2004; Diemand et al. 2004; 



Gao et al.ll2004 iNagai & Kravtsovi i2005i; iMaccio et alji2006i) . 



Studies that include baryons, star formation and cooling, how- 
ever, show that the distribution of galaxies associated with sub- 
halos has a steeper inner profile than the subhalo distribution , 
both at cluster and at galaxy scales dNagai & Kravtsovi 120051: 
iMaccio et al.ll2006l) . For samples of subhalos selected by tidally- 
truncated mass, objects near the halo center lose a greater per- 
centage of their dark matter mass than objects near the virial 
radius. Stellar mass selected samples of satellite galaxies are 
resistant to this effect since baryonic components are located 
in the centers of dark matter subhalos and are tightly bound. 
The observed radial profile of satellites in galaxy-sized halos 
is genera lly more centrally concentrated than the subhalo dis- 
tribution JChen et al.ll2006h but consistent with the dark matter 



profile ("van den Bosch et aljllool IChen et alll2006l) . although 
[Sales & Lambas (2005) have contrary results. 

The color dependence in the spatial distribution of satellite 
galaxies has been studied extensively in angular distributions. 
The angular distribution of satellites is anisotropic and aligned 
with the major axis for red host galaxies and is consistent with 
isotropic for blue hosts (Yang et al. 2006; Agustsson & Brainerg 
i2008l: lAzzaro et al. 2007). These results are consistent with the 
picture where the orientation of the major axis of elliptical galax- 
ies is determined by the direction along which subhalos are 
falling into the galaxy - along filaments. The orientation of spiral 
galaxies is determined by the angular momentum vector, which, 
in simulations, has shown only poor alignment with the mi- 
nor axis of the dark matter halo. A secondary result found in 
these studies is that red satellites of red hosts show a stronger 
anisotropy than blue satellites. This result might be explained 
by a scenario where satellite color is determined by accretion 
time; red satellites were accreted earlier, while blue satellites 
represent recent infall. If the major axis of the host galaxy is 
set early and the orientation of infalling subhalos with respect to 
the host galaxy changes over time, then objects which accreted 
early would show a stronger alignment with respect to the host 
galaxy. 

These observational results have been analyzed by com- 
parison to mock galaxy catalogs derived from semi-analytic 
models - simulating galaxies using the mass accretion histo- 
ries of dark matter halos in a simulation and a variety of as- 
sumption about the physics o f galaxy formation and evolution. 
lAgustsson & BrainerdI (12008 !) use mock galaxy catalogs to sug- 
gest that the difference in degree of anisotropy between blue and 
red host galaxies is real and not due to_differences in the inter- 
loper populations, while iKang et al] (12007 ) come to an oppo- 
site conclusion using a very different technique - t he color de- 
pendence is due to interlopers in the group catalog. IKang et al] 
(2007) also conclude that the difference in alignment between 
red and blue satellites is due to the masses of the satellites; red 
satellites are larger and associated with subhalos which were 
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more massive at the epoch of accretion and which in simula- 
tions are accreted more preferentially along the major axis of 
the halo (Libeskin d et al. 2005; Wang et al. 200 5). The semi- 
analytic model of lAgustsson & BrainerdI (l2008l) suggests that 
the amount of interloper contamination is much greater for blue 
satellites than red satellites (~ 50% to ~ 15%, defining inter- 
lopers as objects with a physical separation greater than 500 kpc 
from the parent galaxy), although the color dependence in the 
anisotropy of satellites is still observed after subtracting inter- 
lopers. 

The color dependen ce in the radial d istribution of satel- 
lites has been studied by S ales et al] (12007) using semi-analytic 
galaxy catalogs constructed from the Millennium Simulation. 
They note that the radial distribution of red satellites is signif- 
icantly more centrally concentrated than the distribution of blue 
satellites. This trend is attributable to accretion time in a scenario 
similar to one explanation for greater anisotropy in the angular 
distribution of red satellites compared to blue satellites: early ac- 
creting satelUtes are stripped of hot gas, stop forming stars, and 
become redder. Early accreters are also likely to orbit closer to 
the center. 

Observationally, the effects of color selection on the radial 
dis tribution of satell i te gala xies has been previously examined 
by ISales & LambasI (|2004|) using data from the Two Degree 
Field Galaxy Redshift Survey (2dFGRS). They find a steeper 
outer slope for satellites of blue parent galaxies than for red pri- 
maries. In addition, red satellites have a distribution that requires 
a larger core radius than blue satellites and they attribute this 
result to correlations between the primary and satellite prop- 
erties. However, in observational searches for satellites, candi- 
date satellites are chosen based upon their projected separation 
and velocity difference from the parent galaxy. Samples of ob- 
jects chosen in this manner are heavily co ntaminated by ob - 
jects which are not satellites - interlopers. In lChen et al.l (|2006|) . 
we discussed the importance of interlopers in making relatively 
unbiased estimates of the projected radial distribution of satel- 
lite galaxies and developed a reliable method of interloper sub- 
traction. Other methods t o account for interlopers h ave been 
investigated. For example, Ivan den Bosch et aLl ( |2004|) exclude 
interlopers us ing an iterative, adaptive selection criterion for 
satellites, and IChenI (l2008l) model satellites and interlopers to- 
gether using a halo occupation distribution (HOD) based ana- 
lytic rn odel for the galaxy correlation function. The IChen et all 
(|2006|) approach has the advantage that it may easily be com- 
pared to previous results. 

In this paper, we add color selection to t he interloper estima- 
tion method described in IChen et al] (|2006|) and apply it to data 
from the Sloan Digital Sky Survey (SDSS) spectroscopic sam- 
ple. We show that samples of red and blue candidate satellites 
have different levels of contamination by interlopers, while the 
model interloper samples for red and blue primaries are similar 
in color distribution. We discuss these results in terms of the en- 
vironmental dependence of satellite galaxies and wider applica- 
tions to the angular distribution of satellite galaxies and to galaxy 
formation models. 



2. Observational Data 

2.1. The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS) includes imaging of 
the the northern Galactic cap in five bands, u,g,r,i,z, down 
to r ~ 22.5 using a dedicated 2.5m telescope at Apache 
Point Observatory in New Mexico in addition to spectro- 



scopic obser vations for a sub s ample of objects from the imag- 
ing catalog (Yorket al."2000'; Fukugita et al. 1996; 'Gunnet alj 
1 998 ; Hoggetal. 2001; Smith etal. 2002; Strauss et al. 200j 
Bla nton et al.ll2003t iGunn eTaH 120061 ITucker et aUMMi . The 
SDSS spectroscopy is carried out using optical fibers posi- 
tioned in pre-drilled holes on a circular plate, with mini- 
mum separation between fibers of 55", the fiber collision dis- 
tance. Reobservations of a field can result in observed spec- 
tra with separations less than the fiber collision distance, down 
to the fiber diameter of 3". The spectroscopic targets are se- 
lected with r-band Petrosian magnitudes r < 17.77 and r- 
band Petrosian half-light surface brightnesses /igo < 24.5 mag 
arcsec"^. An automated pipe line measures the redshifts and 
classifies the redu ced sp ectra dStoughton et alj|2002l: iPier et alj 
2003yivezic et alj|2004l, D. J. Schlegel et al. 2008, in prepara- 



tion)C 

We use the spectrosc opic Main galaxy catalog availa ble as 
Data Release Six (DR6; lAdelman-McCarthv et a"Lll2008l) . cov- 
ering an area of 7425 deg^. Because the SDSS spectroscopy 
is taken through circular plates with a finite number of fibers 
of finite angular size, the spectroscopic completeness varies 
across the survey area. The resulting spectroscopic mask is 
represent ed by a com bination of disks and spherical polygons 
(Tegmar ket al.l |2004|) . Each polygon also contains the com- 
pleteness, a number between and 1 based on the fraction of 
targeted galaxies in that region which were observed. We ap- 
ply this mask to the spectroscopy and include only galaxies 
from regions where completeness is at least 90%. We use r- 
band magnitude s in DR6, built from the NYU Value-Added 
Galaxy Catalog (iBlanton et alJl2005h . normalized to h-1, such 
that Mr - Mo.i^ - 51ogio/z, where Mq.i,. is the absolute magni- 
tude K-corrected to z^O.l (kcorrect v4 . 1 . 4) as described in 
IBlanton & RoweisI (l2007l) . 



2.2. Volume-Limited Galaxy Samples 



Following the procedure in lChen et al.l (l2006l) . we use a volume- 
limited galaxy sample with a depth of 13,500 km s ', corre- 
sponding to the limiting redshift of z = 0.045. This limit is cho- 
sen as a trade-off between the volume of the sample and the ab- 
solute magnitude limit for our satellites, which would need to be 
decreased to brighter magnitudes for more distant primaries. In 
addition, the limiting magnitude sets a minimum separation at 
which fiber collisions become important, which increases with 
distance. To include more distant primaries we would have to 
sacrifice the ability to probe density distributions at small sepa- 
rations. 

From our volume-limited sample, we construct a primary 
sample of isolated host galaxies and a sample of potential satel- 
lites that are projected close to primaries and refer to these two 
samples as the primary sample and the satellite sample. Isolated 
host galaxies are chosen in order to reduce the number of galaxy 
groups selected and eliminate the contamination from satellites 
of galaxy group mem bers. Parameters for the criteria follow 
iPrada et alj ( l2003l) and lChen et al] (120061) and are Hsted in Table 
[Tj The isolation criterion requires that a primary have only neigh- 
bors at least two magnitudes fainter within AR = 0.5 /i ' Mpc 



' We use the reductions of the SDSS spectroscopic data performed 
by D. J. Schlegel et al. (2008, in preparation). The redshifts found are 
identical to the redshifts found by an alternative pipeline used for the 
SDSS Archive Servers (M. SubbaRao et al. 2008, in preparation) over 
99% of the time for Main galaxy sample targets. 
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Table 1. Selection & Isolation Criteria 



Parameters 



Value 



AM, < 2 


0.5 


1000 


AM, > 2 


0.0329 


0.5 


500 


13,500 km s 


1602 


690 


-17.77 



Constraints on primaries M,. < -20 

Isolation criteria: 

Magnitude difference 

Minimum projected distance, hR (A"' Mpc) 

Minimum velocity separation, AV (km s"') 
Satellite sample criteria: 

Magnitude difference from the primary 

Minimum projected distance, 6r Qr^ Mpc) 

Maximum projected distance, Sr (h^' Mpc) 

Maximum velocity separation, Sv (km s"') 
Maximum depth of sample 
Number of isolated primaries 
Number in satellite sample 
Limiting magnitude M, 
Nearby points mock primary criteria: 

Distance from primary Ai^cMi (/;"' Mpc) 1.0-2.0 



and Ay = 1000 km s' Potential satellites of any isolated pri- 
mary must be at least 2 magnitudes fainter than the primary 
and within Sr = 0.5 /i"' Mpc and 6v - 500 km s"'. The max- 
imum absolute magnitude for satellites is set by Mi-n^ - Slogh - 
11.11 - DM - /To.i in the r-band, where the 17.77 is the flux 
limit in this band, DM is the distance modulus, and Kq\ is the 
K-correction at z=0. 1 . This gives a limiting absolute magnitude 
of Mr - -11.11. The satellites are thus limited to the bright- 
est satellite galaxies, ~ O.IL*. In order to avoid biasing from 
a deficit of close pairs of objects, the minimum separation be- 
tween fibers is 32.9 /;"' kpc - the fiber collision separation at the 
redshift of the furthest point in our sample. Finally, we choose 
galaxies that are in areas that are at least 90% complete. For the 
range -23 < Mr < -20, there are 1602 primary galaxies and 
690 objects in the satellite sample with projected radii greater 
than the minimum separation and less than 0.5 /i ' Mpc. 

2.3. Interloper Subtraction 

There is a fraction of objects in our satellite samples that are 
not gravitationally bound to the primaries but are included in 
the sample because of projection effects. Throughout this pa- 
per, we call such objects interlopers (in turn, satellites samples 
without interloper contamination are referred to as true satellite 
samples). I r iterlop er subtraction is discussed in greater detail in 
IChenetalJ (120061) . 

In semi-analytic galaxy catalogs and in dark matter simula- 
tions, interlopers are si gnificant in samples of satellites or DM 
subhalos. For example lAgustsson & BrainerdI (200 8) find in a 
semi-analytic galaxy catalog constructed from the Millennium 
Simulation - using a stricter isolation criteria than we use - that 
the interloper contamination fraction is ~ 30%. In DM-only sim- 
ulations, the fraction of interlopers as a function of projected ra- 
dius rises from a few perce nt at i^ ~50/i~'kp c to nearly 100% at 
R =0.5 /j-^Mpc (see F ig. l. lChenet al.ll2006h . 

IChen et al] ( l2006l) developed and tested several methods of 
subtracting interlopers from the satellite sample statistically. The 
projected surface density of candidate satellites, (2(7?))sat is es- 
timated in bins and normalized by the total number of primary 
galaxies in the sample. Simple methods of interloper subtraction 
estimate the corresponding projected surface density in interlop- 



ers and subtract it from the interloper contaminated surface den- 
sity: <2(/?))i„,,ub. = (2(/?)),,, - <E(/?))i„,. 

A simple and common method of interloper subtraction as- 
sumes a uniform distribution of interlopers in space and esti- 
mates the interloper surface density by counting the number of 
objects that satisfy the satellite criteria but are located around 
empty points in the sky, or mock primary galaxy positions. The 
mock satellite samples created from these points constitute sam- 
ples of pure interlopers: there are no satellites in these samples. 
Hereafter, we generically refer to these estimates of purely inter- 
loper - mock satellite - samples as model interlopers. We then 
fit the surface density profile by a power-law. 



1(R) = AR", 



(1) 



^ For this reason, we only search for primaries within the subset of 
velocities 1000 to 12,500 km s"' . 



with slope, a, and normalization, A. In DM simulations, the dif- 
ference in slope between the satellite sample with no interloper 
subtraction and the true satellite sample is Aa ~ 0.5. The differ- 
ence between the satellite sample and the interloper-subtracted 
sample using this simple estimate is Aa ~ 0.1, a vast underes- 
timate of the fraction of interlopers. This simple interloper sub- 
traction method fails because it oversamples voids compared to 
clustered areas - there are a lot more isolated, empty points in 
voids than elsewhere. Most galaxies, however, reside in clustered 
areas, even for our sample of isolated galaxies. 

To account for this clustering effect, our preferred method - 
the 'nearby points' method - samples the environments of our 
primary galaxies. We pick mock primary galaxy positions that 
are between 1 /;"' Mpc and 2 /i"' Mpc from a real primary: the 
minimum separation avoids sampling real satellites. All mock 
primary points satisfy the same isolation criteria as our sample 
of primaries, and the size of this mock primary sample is 20 
times the number of primary galaxiesQ 

Model interlopers from the area near actual primaries better 
approximate the interloper contamination in the satellite sample. 
For example, mock galaxy catalogs show that the sample of true 
interlopers contains objects with velocities correlated to the pri- 
mary: the distribution of relative velocities between interloper 
and primary is not constant (see Fig. 3. IChen et alJ i2006). This 
structure in relative velocities of true interlopers is reproduced 
in the relative velocities of model interlopers with respect to pri- 
maries using the nearby points method and cannot be reproduced 
using the simple interloper method outlined above. In previous 
tests on mock galaxy catalogs, the nearby points method returns 
a best-fit po wer-law slope with in ~ 0. 1 of the true satellite dis- 
tribution. In IChen et aiT ( l2006l) . for a set of SDSS galaxies us- 
ing the same isolation and satellite criteria as used in this paper, 
the nearby points method returns a best-fit power-law slope of 
Q' = -1.58+0.11. The eiTors quoted are statistical and do not in- 
clude the Aa ~ 0.1 bias inherent to the method. Accounting for 
the systematic error would give cc w -1.7 + 0.1. Using our larger 
DR6 data set, we calculate the surface density of objects in six 
bins of width 0.067 /;"' Mpc, and we find a consistent result for 
the nearby points method, a = -1.52 + 0.07 (see T able[2l). These 
results are also consistent with 'S ales et al.l (12007 ') who estimate 
the radial distribution of satellite galaxies using semi-analytic 
galaxy catalog constructed from the Millennium Simulation and 
finda^ 1.55 ±0.08. 

As a final note, we might expect that the interloper- 
subtracte d su r face d ensity profile would be best fit by the 
Navarro et al.l d 19971 NFW) profile that is used to describe the 
distribution of matter in dark matter halos. However, our mini- 
mum separation is on the order of the expected scale radius for 



In our data there is a total of 3325 model interlopers. 
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Fig. 1. The color distribution for the primary sample (solid), 
satellite sample (dotted), and model interloper sample (dashed). 
Red satellites are defined to have Mg - Mr > 0.65, while red 
primaries have colors Mg - M,- > 0.9. Poisson errors are shown 
for the primary and satellite sample. The corresponding Poisson 
errors for the model interloper sample are significantly smaller 
than for the satellite sample. 



our parent galaxy halos, making fitting for the scale radius im- 
practical. 

3. Results 

The Mg - Mr color distributions for the primary, satellite, and 
model interloper samples are shown in Fig.[T] The primary sam- 
ple is redder than the satellite sample with a large red peak and a 
bluer tail, while the satellite sample has a vaguely bimodal distri- 
bution, split at Mg-Mr ~ 0.65. The model interloper distribution 
appears tilted to bluer objects compared with the distribution of 
the satellite sample. 

3.1. Red and Blue Satellites 

The differences in color distribution between the satellite sample 
and the model interloper sample suggests that a greater fraction 
of blue objects in the satellite sample will be interlopers than of 
red objects. Defining red satellites to Mg-Mr > 0.65, of the 690 
objects in the satellite sample, 285 are red and 405 are blue. By 
comparison to the model interloper sample, the fraction of the 
objects that are interlopers in the whole sample is 24%, while 
the percentages for the red and blue satellite sample are 16% 
and 30%, respectively. 

This difference affects the estimates of the slope of the 
power-law fit to the radial distribution. In Figure|2] the biasing of 
the best-fit power-law is shown for all satellites and blue and red 
satellites. The distribution of blue satellites is significantly flatter 
than the distribution of red satellites, and the distribution of blue 
satellites is more effected by interlopers than the distribution of 
red satellites. 

In the interloper-contaminated satellite samples, the best-fit 
slopes for blue and red satellites (shown in Table 2 and Figure[3j 
are inconsistent, with a much shallower slope for blue satellites. 
After interloper subtraction, although the difference in slopes is 
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Fig. 2. The surface density of the satellite sample (squares) and 
the interloper-subtracted satellite sample (circles) for all satel- 
lites (top), blue satellites (center), and red satellites (bottom). 
The best fit power-laws are also plotted with solid lines for the 
satellite sample and dotted lines for the interloper-subtracted 
sample. 



smaller, the blue slope is still significantly shallower than the red 
slope, and the slopes are inconsistent given both the marginal- 
ized errors and their Icr confidence intervals (see Fig. [3). The 
projected radial distribution of red satellites is as steep as might 
be expected for the dark matter density distribution in halos 
which host galaxies such as are found in our primary sample. 
The distribution of blue satellites, on the other hand, is as shal- 
low as might be expected for the subhalo distribut ion in the same 
host halos (see, for comparison. lChen et al.ll20()6l) . 
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Table 2. Best-Fit Power-Law Slopes to the Surface Density 
Profile of Satellites 



Input Data 



Satellite Sample 



Interloper Subtracted 



all 


-1.24 + 0.06 


-1.52 + 0.07 


blue satellites 
bright 
faint 


-1.01+0.08 


-1.32 + 0.10 
-1.42 + 0.15 
-1.25 + 0.13 


red satellites 
bright 
faint 


-1.59 ±0.09 


-1.80 + 0.10 
-1.78 + 0.14 
-1.86 + 0.15 


blue primaries 
bright 
faint 


-1.18 + 0.10 


-1.54 + 0.12 
-1.67 + 0.15 
-1.35 + 0.20 


red primaries 
bright 
faint 


_i 97+0.07 


-1.50 + 0.09 

-1.34 + 0.11 

-1 82+"'5 

^•°^-0.14 



-1.5 



-2 



■w 




0.02 0.04 0.06 0.08 
A 



0.1 



Fig. 3. The 68% confidence intervals for the normalization (A) 
and slope (a) of the power-law, I.{R) = AR". Shaded contours 
are for blue satellites; solid line contours are for red satellites. 
The contour for interloper-subtracted samples are in each case at 
steeper slopes than their satellite sample counterparts. 



In general, we expect blue satellites to be fainter than red 
satellites, so it is useful to attempt to disentangle the effects of 
luminosity from those of color. While Chen et al. (2006) tests 
luminosity dependence (which is further discussed in Section 
m, they do not have suflicient numbers of objects to test color 
and luminosity together. Figure |4] shows the luminosity distri- 
butions for the blue and red satellites and primaries. We split 
each color sample of satellites into a faint and a bright sample 
at Mr = -18.3. Table |2] shows that the interloper-subtracted sur- 
face density profiles of faint and bright samples are consistent 
for color-selected samples, i.e., color is the dominant attribute in 
the radial distribution of satellites. 



3.2. Red and Blue Primaries 

We split the 1602 primary galaxies by color at Mg - M,. - 0.9. 
The 591 red primary galaxies have 358 objects in their satellite 
sample, while the 1011 blue primary galaxies have 332 candi- 
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Fig. 4. The luminosity distribution with corresponding Poisson 
errors of primaries {top) and objects in the satellite sample {bot- 
tom), where red objects and blue objects are shown separately, 
in dotted and solid lines respectively. Bright and faint satellites 
are split at M,- = -18.3, while bright and faint primaries are split 
at Mr = -20.8. 



date satellites. The interloper percentages for the red and blue 
primary sample are 21% and 28%, respectively. Red primaries 
of Mg - Mr > 0.9 have a larger fraction of red objects in their 
satellite sample than than blue primaries (52% to 30%) and, sub- 
sequently, more true satellites (see Fig.|5]). This can also be seen 
in Fig.|6l where the amplitude of the satellite profile of red pri- 
maries is greater than that of the blue primary profile. After inter- 
loper subtraction, the slopes of the density profiles of satellites 
for blue and red primaries are consistent. 

Figure |4] additionally shows that red primaries are in gen- 
eral brighter than blue primaries. Splitting the color-selected pri- 
maries into faint and bright samples at M, = -20.8, Table |2] 
shows significant differences for the distribution of satellites of 
faint and bright primaries, depending on the color of the primary. 
Red, brighter primaries have an interloper-subtracted satellite 
profile that is shallower than the distribution of satellites in red, 
fainter primaries, a = -1.34 + 0.11 to a - -1.82^^ J^. This is 
not unexpected; from numerical simulations we expect the satel- 
lite distribution to scale with the mass distribution of the primary. 
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Fig. 5. The color distribution for the satellite sample and model 
interloper sample for red primaries (thick and thin solid lines, 
respectively) and the color distribution for the satellite sample 
and model interloper sample for blue primaries (thick and thin 
dashed lines, respectively). 



and larger primaries have smaller concentrations leading to shal- 
lower slopes at the radii at which we measure. Interestingly, 
blue primaries show the opposite dependence: the distribution 
of satellites around brighter primaries seems to have a steeper 
slope than around fainter primaries. This trend is accompanied 
by a change in the fraction of red satellites: 38% of the satel- 
lite sample for bright, blue primaries are red, while 23% of the 
corresponding sample for faint, blue primaries are red. The er- 
rors, however, are nearly as large as the discrepancy between the 
samples, and larger samples will be required to confirm this re- 
sult. In addition, the slope for satellites of bright, blue primaries 
is steeper than that for bright, red primaries, which may reflect 
color dependence in using luminosity as a proxy for mass. 

Similar results can be seen if we split both the primary and 
the satellite samples by color In Table[3] we show that red satel- 
lites around blue primaries have a larger interloper fraction and 
a steeper power-law slope than found for red satellites of red 
primaries. The difference in power-law slope can be attributed 
to the same luminosity dependence discussed previously: red 
primaries are more massive and have mass and satellite dis- 
tributions described by smaller concentrations. Interloper con- 
tamination increases from small separations to large, and for 
larger primaries, the range of radii we probe preferentially sam- 
ple areas with smaller level of interloper contamination. On the 
other hand, for blue satellites of red and blue primaries, we find 
very similar interloper fractions and shallow power-law slopes. 
Unfortunately, in all cases, statistics are poor and future larger 
samples will be required for statistical significance to be as- 
cribed. 



3.3. Environmental Dependence 

Despite the difference in satellite samples, the color distribution 
of the model interloper samples of red and blue primaries are 
similar to each other and to the color distribution of the satel- 
lite sample of blue primaries (see Fig.|5]l. Restated, although red 
primaries are found in more clustered environments, the envi- 
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Fig. 6. The surface density of the satellite sample (squares) and 
the interloper-subtracted satellite sample (circles) for satellites 
of blue primaries (top) and of red primaries (bottom). The best 
fit power-laws are also plotted with solid lines for the satellite 
sample and dotted lines for the interloper-subtracted sample. 

Table 3. Primaries and Satellites Selected by Color 



Input Data 


Number in 


Interloper 


Power-Law 


(Primary-Satellite) 


Satellite Sample 


Fraction 


Slope 


blue - blue 


232 


0.30 


-1.33 ±0.15 


blue - red 


100 


0.22 


-1 98+01'' 


red - blue 


173 


0.30 


-1.26 ±0.14 


red - red 


185 


0.12 


-1.71 ±0.12 



ronments of red and blue primaries are not noticeably different 
as measured by the color distribution of model interlopers. 

Red primaries live in more clustered environments than blue 
primaries as measured in the average surface number density of 
model interlopers: 0.159 h^ Mpc"^ for red primaries and 0.117 
h^ Mpc"^ for blue primaries. In the previous section, we noted 
that these slightly different environments have similar color dis- 
tributions of faint objects, even though we do not use fixed lumi- 
nosity criteria for these objects. We test one less clustered envi- 
ronment, sampling isolated points that lie outside of a 2 /i"' Mpc 
radius from our isolated galaxies (outside the criteria for our in- 
terloper subtraction method). Here, the average surface density 
is 0.042 h^ Mpc"^. All three model interloper distributions are 
plotted in Fig. |7] in which all environments have similar color 
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Fig. 7. The color distribution for the model interloper samples 
for red primaries (solid), blue primaries (dashed), and for sam- 
pling isolated points greater than 2/1"' Mpc from primaries (dot- 
ted). The average surface number density of objects in these 
samples is 0.159, 0.117, and 0.0415 h^ Mpc I 



dependences with a possible trend to bluer objects in less dense 
fields. 



4. Conclusions 



In IChen et al.l (l2006l) . we constrained the projected radial dis- 
tribution for isolated galaxies and found that their power-law 
slopes are steeper than the expected slopes for the distribution 
of dark matter subhalos and may be as steep as the density pro- 
file of the host dark matter halos. We reproduce this result with 
a survey area that is ~ 50% larger. However, the distribution of 
satellite galaxies shows some significant dependence on color. 

When samples of candidate satellites are split by color, we 
see that blue objects are more likely to be interlopers than red 
objects. The observed estimated interloper contamination of red 
objects in the satellite sam ple is 16%, whi le that of blue ob- 
jects is 30%. Agustsson & Brainerd I (l 2008 h produce a differ- 
ence of 15% and 50%, respectively, in the interloper contami- 
nation of red and blue satellites samples, using a semi-analytic 
galaxy formation model. It is, then, not appropriate to assume 
that red and blue satellites have the same level of interloper con- 
tamination and this must be taken into account in testing the 
color dependence in the angular distribution of satellite galaxies. 
This result also suggests that the prevalence of faint red galax- 
ies could be developed into a method to find small groups of 
galaxies in a manner similar to the method by which the red se- 
quence of early-type g alaxies is used to find galaxy clusters (e.g., 
iGladders & YeSwOti . 

Both with and without interloper subtraction, the radial pro- 
file of blue satellites is significantly shallower than that of red 
satellites. The best-fit power-law slope of interloper-subtracted 
blue satelHtes is a = 1.34 + 0.10 and the that of red satellites is 
a = -1.80 + 0.10. T his result is consistent with the trend found 
bv lSales et al.l (1200 7*) who test color dependence in semi-analytic 
galaxy catalogs constructed from the Millennium Simulation, 
finding a more centrally concentrated radial distribution of red 
satellites than of blue satellites. This trend is repeated when they 



select by a proxy for accretion time - whether the satellite retains 
a DM halo or whether it has been tidally destroyed. 

Blue satellite s are gener a lly fa inter than red satellites. 
Correspondingly, IChen et aTl (l2006h found that the best-fit 
power-law slope for bright satellites is steeper than that of faint 
satellites (cut at Mr = -18.28) - although without statistical sig- 
nificance - in a volume-limited sample. In a flux-limited sample, 
the reverse relation was found; however, this sample used a lu- 
minosity cut as faint as the faintest satellites in our sample, at 
Mr - -17.76, and so is not directly comparable. When we cut 
the satellites samples by color and luminosity, we find that the 
dominant effect is from color; the best-fit power-law slopes of 
bright and faint samples of red satellites are consistent as are 
those of bright and faint blue satellites. 

The power-law slope for red satellites is as steep as might be 
expected for the dark matter density distribution of halos which 
host galaxies like those found in our primary sample. On the 
other hand, the power-law slope for blue satellites is as shal- 
low as the expected subhalo distribution. The shallowness of the 
subhalo profile is attributed to tidal stripping. This is unlikely to 
effect satellite galaxies, since they are located at the centers of 
the dark matter subhalos. The shallower profile for blue satel- 
lites as compared to red satellites, then, might be interpreted as 
consistent with the scenario where satellite color is determined 
by accretion time: red satellites were accreted earlier, while blue 
satellites represent more recent infall. Satellites are expected to 
be redder in the inner regions of parent halos due to environ- 
mental processes that shut off star formation (ram pressure strip- 
ping, strangulation, etc.). This morphological segregation has 
also been observed in more massive clusters and in the fainter 
satellites found in the Local Group. 

When splitting the satellite sample by primary color, red 
primaries have a significantly larger fraction of red satellites 
and somewhat smaller interloper fractions than blue primaries. 
After interloper subtraction, the best-fit power-law slopes of 
satellites of red and blue primaries are consistent within errors. 
Red primaries are, on average, more lum inous than blue galax- 
ies. Correspondingly, ' Chen et all ( 1200 6') found that the best-fit 
power-law slope for satelUtes bright primaries is consistent with 
that of satellites of faint primaries (cut at Mr = -21). When we 
cut the primary sample by luminosity and color, the trend with 
luminosity is different for red and blue primaries. The slope of 
satellites of bright, red primaries is shallower than that of faint, 
red primaries, a relation probably dominated by the mass of the 
parent halo as brighter primary galaxies reside in bigger par- 
ent halos which have mass distributions characterized by smaller 
concentrations. The slope of satellites of bright, blue primaries 
is steeper than that of faint, blue primaries, as the fraction of red 
satellites drops with primary luminosity. 

The color distribution of objects in model interloper sam- 
ples are similar, regardless of their environment (as measured 
by average surface density). While robust conclusions cannot 
be drawn as to how similar these distributions are, it suggests 
that there are fundamental differences between satellite galax- 
ies and faint galaxies in the field. Intriguingly, this blue-tilted 
color distribution also resembles that of the satellites sample of 
blue primaries. Better understanding of the processes that effect 
the color of faint objects and the radial distribution of satellite 
galaxies will require further studies. 
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